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Many models for physics beyond the Standard Model predict lepton-flavour violating decays of charged leptons 
at a level which may become observable very soon. We investigate the decays of a r-lepton into three charged 
leptons (r ^ IW , 6''^ = e, ^i) in a generic way. Using effective-field-theory methods, the relevant operators are 
classified according to their chirality structure. For each case, we work out the Dalitz plots for the decay distribu- 
tions, including interference terms which arise from four-lepton operators and radiatively induced processes. We 
discuss phenomenological implications, in particular the potential to distinguish different new physics models. 



1. Introduction 

Over the last few years lepton flavour viola- 
tion (LFV) has become a hot topic. On the one 
hand, neutrino oscillations have been discovered 
recently, which is already a clear signal for physics 
beyond the Standard Model (SM). Incorporating 
these effects into a minimal extension of the SM, 
LFV for the charged leptons is predicted, but 
however at a completely unobservable level. On 
the other hand, many models forcast measurable 
new physics (NP) effects at the TeV scale, which 
can be investigated soon with the upcoming ex- 
periments. These extensions of the SM predict 
LFV at much higher rates, which, in some cases, 
may already be in conflict with existing experi- 
mental bounds [1] (see also for a recent sum- 
mary of i?-factory results). With the advent of 
new experimental facilities |3] (see also U) the 
current bounds will be pushed further, if not a 
discovery will be made. 

To identify effects beyond the SM one needs, 
among others, reliable bounds on LFV decay 
modes. This requires knowledge on the phase 
space dependence. Furthermore if a discovery 
is made, one can directly compare different NP 
models. But due to the huge number of different 
models, this is quite cumbersome. On the other 
hand a model-independent approach allows us to 
classify the effects in a general way. This results 
can then in principle be used to compare or re- 
strict different models. 

All this models which predict LFV r decays of 



the form r ^ U'i" with iJ'J" ^ e, fi |5ieil7l8l9l 
110111112113114115116] will eventually match onto a 
set of local four-fermion operators or radiative op- 
erators, the latter mediating r £7* with subse- 
quent decay of the (virtual) photon into a charged 
lepton pair. 

We persue a bottom-up approach, where we re- 
gard the SM as an effective-field-theory (EFT). 
The NP effects are expanded in terms of 1/A, 
where A is the scale corresponding to the NP 
scenario. This allows us to consider all possi- 
ble four-fermion and radiative operators with ar- 
bitrary coupling constants, which can be deter- 
mined by studying the decay distributions of the 
three leptons in the final state. NP models are 
then distinguished by using the fact that in dif- 
ferent models different operators dominate. Even 
if no signal events are found, such a study of the 
decay distributions is necessary to determine the 
efficiency of an experiment and hence to extract 
reliable limits. 

At the LHC experiments it will be possible to 
detect LFV decays of a r lepton, especially into 
channels with three leptons. The signal r — > 3/z 
will be one of the cleanest signatures [17 ; we will 
therefore concentrate on this in the following. 

First we classify the general set of effective op- 
erators and formulate the corresponding interac- 
tion theory. After pointing out the kinematical 
variables, we show the results for the decay dis- 
tributions. In the end wc shortly comment on the 
unknown couplings in minimal ffavour violating 
scenarios and present a conclusion. 
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2. Theoretical Basics 

2.1. Operator Analysis 

We construct the EFT by assuming that some 
contribution from unknown physics beyond the 
SM at a high scale A induces LFV processes]^ 
At the electro-weak scale these LFV interactions 
manifest themselves in higher dimensional op- 
erators, which have to be compatible with the 
SU{2)l X U{1)y gauge symmetry of the SM. The 
leading operators for processes involving charged 
leptons will be of dimension 6j^ We will construct 
these operators in the following. 

To respect the SM gauge symmetry we group 
the left-handed leptons in a SU{2)l doublet, 
while the right-handed charged leptons (which are 
singlets under SU(2)l) are put into an incom- 
plete doublet, as a reminiscent of a right-handed 
SU{2)ji related to custodial symmetry. Writing 
also the Higgs boson in matrix form, we have 



L = 





K L ) 


) , R = 





H 



1 / + /iQ + »Xo 



72 



). (1) 



For simplicity we have suppressed the family in- 
dices, which will be specified once we consider a 
particular decay mode. In terms of the fields de- 
fined in ([l]), the list of operators relevant for the 
lepton flavour violating t W decays reads 

dimension 6 leptonic: 

01 = {Lj,,L)iL^''L) (2a) 

02 = {LT''^^L){LT'^rL) (2b) 

03 = {Ri^R){RrR) (2c) 

04 = {RiMLl'^L) (2d) 

dimension 6 radiative: 

i?i = g'{LHa^,R)B^''' (2e) 
i?2 = giLT''Ha^,R)W''''^'' (2f) 

^Notice that, in general, the scale associated to lepton- 
flavour violation is independent of the scale related to 
lepton- number violation, Aln- 

^Dimension 8 operators, which involve scalar and tensor 
currents, are further suppressed by j . The suppres- 
sion in 2HDM could be reduced due to large tan /3 effects. 



where is the J7(l)y gauge field, are the 
SU{2)l gauge fields and g and g' are the cor- 
responding gauge couplings. In this analysis we 
concentrate only on operators which contribute 
at tree-level; operators which are bi-linear in the 
lepton fields and contribute at the loop-level, can 
be found e.g. in j 7|9|18j . The most general effec- 
tive Hamiltonian at the electro-weak scale is then 
obtained by summing over these operators, mul- 
tiplied by arbitrary coefficients for every ffavour 
combination. In a particular new physics sce- 
nario, these coefficients should be obtained by 
matching at the new physics scale A and evolving 
down to the scale M\n/ within the SM which is 
considered as an effective theory (ET). 

To describe LFV decays of a r lepton into three 
charged leptons, we have to construct the effective 
interaction at the scale of the r lepton. To this 
end we integrate out the weak gauge bosons and 
the Higgs. In the following we will focus on t~ 
decays; the decay distributions for t+ decays are 
identical. On the level of four-fermion operators 
with dimension 6, we obtain the same structures 
as in ( 2apd I. Projecting on charged leptons only, 
we see that O2 becomes equivalent to Oi, and 
both match onto a purely left-handed operator 



H 



{LL){LL) _ JLLKLL) {(■Llt^TLWLl^'i'l) 



off 



9V 



A2 



(3) 



where here and in what follows the superscript of 
the coupling denotes the combinations of chiral- 
ities involved and the subscript denotes the rele- 
vant Dirac structure. Likewise, the operator O3 
corresponds to a purely right-handed interaction 



H 



(R.R)(B.R) _ JRB.){RR) (^fl7pTi?) (^K7^^fl) 



off 



^9v 



A2 



(4) 



while we get a mixed term from the operator O4 



H. 



{LL){RR) _ {LL){RR) {^Ll t,TL){i'ii"f^ C-'b) 



off 



--9v 



A2 



^ (RR,)iLL) iiR7^Tn)i£'^r£l) 

+9v ■ ^ ' 

In case of the dimension 6 radiative operators 
we are only interested in the neutral current com- 
ponent coupling to a charged lepton pair, which 
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mediates the LFV decay into three charged lep- 
tons. Switching to the physical fields accounting 
for the neutral current component, the photon 
and fields, we have to integrate out both the 
Higgs and the Zq. In doing so we obtain a ra- 
diative operator with a photon, as well as a four- 
fermion contribution from Zq exchange. We find 
that this operator related to Zq exchange is sup- 
pressed relative to the photon exchange and four- 
fermion operators by the small Yukawa coupling 
of the T lepton. Thus, only the photonic contri- 
bution has to be taken into account. For this we 
obtain at the scale m,^ 



rrrad 



47r A2 



h.s 



(L Ft) ( Ft L) 

where g^^^ and g^^^ denote the two possible 
chirality combinations]^ The matrix element for 
r £i'i' becomes 



{££'£' \H^ 



rad I 
off I 



V q 



A2 g2 



{fe£'\ {h{-^a,.)r,) [t'lH') 



(7) 



where q is the momentum transfer through the 
photon. This momentum transfer is proportional 
to the lepton masses, and thus this contribution 
scales as l/iyA"^) where y is a Yukawa coupling of 
the leptons, which would lead to an enhancement 
unless an additional Yukawa coupling appears in 
the numerator as, for instance, in minimal flavour 
violation (MFV). 

Four different Dirac structures can occure if 
only the four-lepton operators contribute. Two 
additional Dirac structures are possible for de- 
cays of the form (! where also the 
radiative operator has to be taken into account. 



^ Here, for simplicity, we neglected possible form factor 
effects for decays into virtual photons from long-distance 
lepton or quark loops. In the most general case, the r — > 
£7* vertex could be parametrized as 



47r A2 



where a\"J^'' {Q) 



^ S^ad ' '^nd f!^'J\0) = 0, see e.g. [Tg]. 



The corresponding coupling constants 



(LL)(LL) 

yv ' 



(LL)(R.R) 

yv 1 



(RM)(LL) 

yv 



of the four-lepton operators and additionally 



{LR) 
i/rad ' 



i/rad 



from the radiative operators are matrices in lep- 
ton flavour space. 

2.2. Decay Modes 

There are in total six different decay modes of 
the to consider 



e e e 



(8a 
(8b 
(8c 
(8d 
(8e 

(8f: 



Notice that ( 8a - 8d I contain two identical parti- 
cles (e~e~ or fi^fj,^) in the flnal state for which 
we have to take into account the Pauli principle 
(PP), whereas ([8e| -t- [8f| do not. Moreover, only 
( [8a| |8b[ [8e| [8f| receive contributions from the ra- 
diative operators (|6| via 



£-7* 



in. 



Therefore we calculate the Dalitz distributions 
d^r/(d?TiJ_dm2_). We neglect interference terms 
with different helicities, because of the suppres- 
sion due to small lepton masses, and take into 
account only the leading interference term. 



3. Results 

The Dalitz distributions are visualised in form 
of density plots, made in such a way that the high- 
est value is displayed as black. Here we concen- 
trate on the difference between four-lepton and 
radiative operators and show only the results to 
T~ fi^ fi^ fi^ . Details to all different possibil- 
ities can be found in [2U]. The kinematical vari- 
ables are then given by 



'12 



TO I 



{Pu- +Pu+? 



4 



ml_ I GeV^ 



Figure 1. Dalitz Distribution for 
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Figure 3. Dalitz Distribution for 



> 

O 1.5 



1 1.5 ; 
/ GeV^ 



conibmation is given by d 



ure,-2.j Simulated 100 events for 
e tnird conibmation is given by 



d^r 



amSo dmr 



The calculation of the decay distribution for the 
vector current, where all particles axe left-handed, 
results in 



^2p(LL)(LL) 

dm^g dTO?2 



|^(Lf,L-)(L^L^')|2 



A4 



■,2 \2 



(2m22 



256 TT-^ TO? 



(9) 



The distribution Q for the vector current with all 
particles left-handed is displayed in fig. [l] in form 
of a density plot. Furthermore a simulation of 100 
events based on (|9| can be found in fig. [2j One 
finds a rather fiat distribution over the allowed 
phase space region. For the radiative transition 



> 
1.5 



Figure 4. Simulated 100 events for 
we find the distribution 



^ ^rad ^ 2 



TO 



128 7r3 TO? 



A4 



'13 



2m22 



3to2 



128 7r3 m? 



TO^ (to^ - Smlml + 2mt) 



128 7r3 TO^g; 



^23) 



'12 



'-13 ■ 



'23 



6to2(to2 -f to2) 



256 TT^ ; 



'■23 ' 



(10) 



The corresponding density plot can be found in 
fig. [3] Again we have made a simulation, see fig.|4] 
Due to the photon pole we find an enhancement in 
the corners of the phase-space for the two possible 
combinations (because of PP) of the fJ.~fJ.^ pair. 

Finally we have considered the leading interfer- 
ence term, between the radiative transition, cou- 
pling to a left-handed pair of muons and the four- 
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lepton operator with all particles left-handed. 
The suppression is given by tjIt/v, at least in 
MFV scenarios. The distribution is given in (111 



and plotted in fig. [5] The simulation can be found 
in fig. [6] 



.2y(LL)(LL) 
mix 



'23 



dr7i^2 



2t;Ke[gv ffrad \ 

A4 



'12 



3m2 



ml{ml~ml){rnl^ + rril^) 



64 TT'^ ml 



■7-3 2 m 2 'y^o 2 



128 TT^ ml 3 



(11) 
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Figure 5. Dalitz Distribution for ^^-r"^^'! 



SU{3)l X SU{3)En is described by two spurion 
fields 

A = - diag(TOe, m^.mr) 



A 



9v 



LN 



U* d\a,g{m^^ , ^ ^1^2 ) 



where A ~ (3, 3) describes the SM Yukawa cou- 
plings of the charged leptons, and ~ (6, 1) is 
related to a dim-5 lepton-number violating term 



-Maj 



-1- (N-gN) 



2Aln 

The quantum numbers of 



N : 



1 



(12) 



(13) 



are vanishing under the complete SM gauge 
group. Both spurion fields are only parametrised 
by the lepton masses and the PMNS matrix U. 
To render the LFV operators gauge invariant, we 
have to insert appropriate numbers of these spu- 
rions. Here we are only considering the mini- 
mal insertion. For one right-handed field we need 
an insertion of A, which produces a suppression 
of mi/v. The four-lepton operators U U L\L1 
need at least two insertions of g^. The possible 
flavour structures can be read off from the reduc- 
tion of the S'[/(3)l tensor product for g^, and g\, 



6x6 = 1 



27. 



(14) 



> 
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Figure 6. Simulated 100 events for ^^^i 

3.1. Comment on MFV scenarios 

Starting point is a scenario with the minimal 
field content, as present in the SM. The nec- 
cessary breaking of the lepton flavour symmetry 



The singlet term corresponds to the trace of 
tr[gjgi,] and does not induce flavour transitions 
at all. The octett A = A''' describes flavour tran- 
sitions from 2-lepton as well as from 4-lepton op- 
erators. The 27plet combination G^!j accounts 
only for the 4-lepton operators. The dominat- 
ing flavour coeflicents for the relevant operators 
mediating LFV r decays in MFV are given by 



(LfcL')(L,L^) 
ilV 

ilv 

i/rad 



2ciAf (5j + C2G, 
C3Af4 

c4A'''Ai . 



kl 



(15a) 
(15b) 
(15c) 



Therefore, assuming a MFV scenario, we can re- 
late the couplings to the lepton masses and neu- 



trino parameters via the relations (15). 
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4. Summary and Conclusion 

LFV T decays provides an important test of 
new physics contributions against the standard 
model, since many new physics models allow for a 
huge enhancement of this decay channels. There- 
fore a great potential to falsify the SM in the up- 
coming experiments is given. We have discussed 
LFV T £££ processes in a model independent 
way. For our analysis we regard the SM as an 
effective theory and classify different operators 
within this effective theory. The operators dif- 
fer in their chirality and Dirac structure. We 
have calculated the Dalitz distributions for these 
different operators including the leading interfer- 
ence term. It turns out that four-lepton operators 
show a rather uniform behaviour, whereas the de- 
cay distributions of radiative operators show large 
enhancement in the corners of the phase space. 
Therefore this analysis allows to distinguish dif- 
ferent models with already only a few number of 
events, because different models give rather dif- 
ferent predictions for the relative size of radiative 
and four-fermion operators. This clearly shows 
that using Dalitz plots for such processes is a com- 
plementary tool to the sole analysis of the total 
decay rates. Additionaly the knowledge on the 
phase space distributions allows to improve the 
bounds on this decay modes by taking into ac- 
count the different structures. 
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